Mesoscopic bound on anisotropy in itinerant ferromagnets 
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We calculate the anisotropy energy of a single-domain ferromagnetic particle in which the only 
source of anisotropy is the presence of non-magnetic impurities. We find that such anisotropy takes 
the form of combined easy-axis and easy-plane anisotropies, with random orientations of the axes. 
Typically the anisotropy energy is of order N^^^h/rso, where A*' is the number of electrons in the 
ferromagnetic particle and Tso is the spin-orbit time. 

PACS numbers: 75.75.-|-a, 73.22.-f, 75.30.Gw 



Not all magnetization directions are equal for a single- 
domain ferromagnet. Through magnetic dipole interac- 
tions and spin-orbit coupling, the shape of a ferromag- 
netic particle and the underlying crystal lattice lead to 
preferred directions for the magnetization Q . In the ab- 
sence of these anisotropy sources, how isotropic can a 
ferromagnetic particle be? We'll answer this question for 
an itinerant ferromagnet with non-magnetic impurities. 

For shape or magnetocrystalline anisotropy the energy 
gain for a magnetization pointing along a certain 'easy 
axis' or inside an 'easy plane' is proportional to the num- 
ber of electrons N contributing to the magnetic moment. 
Clearly, for the question asked here, there is no such 
bulk effect. However, for small itinerant ferromagnets, 
non-magnetic disorder or irregularities in the boundary 
of the particle break the rotation symmetry on the mi- 
croscopic scale, thus creating a 'random' anisotropy in 
the orbital part of electron wavefunctions. This orbital 
anisotropy affects the magnetization through spin-orbit 
coupling. The resulting 'mesoscopic' anisotropy does not 
scale proportional to iV, but it is a fundamental lower 
bound on the anisotropy energy of a ferromagnetic parti- 
cle. It can be an important contribution in ferromagnetic 
nanoparticles or in disordered ferromagnetic alloys. 

Recently Gueron et al. and Deshmukh et al. 
studied the excitation spectra of ferromagnetic parti- 
cles in the nm size range. The excitation spectrum 
provides an indirect method to measure the magnetic 
anisotropy energy 4, 5, 6, f^. Although shape domi- 
nates the magnetic anisotropy in these nanoparticles, the 
experiments revealed a significant mesoscopic contribu- 
tion to the anisotropy 0, 13- Mesoscopic contributions 
to the magnetic anisotropy energy have also been con- 
sidered theoretically by Cehovin et al. Q and by Usaj 
and Baranger Motivated by the experiments of Refs. 
0, m, these authors studied how the addition of a sin- 
gle electron changes the magnetic anisotropy energy in 
a ferromagnetic particle with a systematic source of uni- 
axial anisotropy. The question we ask here is different: 
it is about the mesoscopic lower limit of the magnetic 
anisotropy when the magnetic anisotropy energy is dom- 
inated by mesoscopic effects. This may be appropriate for 
the smallest ferromagnetic nanoparticles with low mag- 



netocrystalline or shape anisotropy. 

We consider an ensemble of ferromagnetic particles, 
with equal shape, size and concentration of non-magnetic 
impurities, but with a different configuration of impuri- 
ties for each particle. All particles in the ensemble are 
roughly spherical and have no notable magnetocrystalline 
anisotropy. For such an ensemble, the average free energy 
F does not depend on the direction of the magnetization 
m. Hence, the magnitude of the 'mesoscopic' anisotropy 
must be characterized through its sample-to-sample fluc- 
tuations, which can be inferred from the correlation func- 
tion 



C{e) = (F(mi)F(m2)) 



(1) 



where 6 is the angle between the magnetization directions 
mi and m2. 

In principle, the ferromagnetic exchange field must be 
calculated self consistently and at a fixed number of elec- 
trons, taking into account the disorder. We simplify the 
calculation by fixing the exchange field a priori (parallel 
to m) , and by performing the average at a fixed chemical 
potential /x, rather than a fixed number of particles. The 
latter simplification does not si gnif icantly affect the fluc- 
tuations of the free energy [ifl, llll ll^l • The system un- 
der consideration is then described by the single-electron 
Hamiltonian 

n^Ho + V, no^^~Ezm-(T-^i, (2) 



2m 

where Ez — fiBBcx is the Zeeman energy corresponding 
the exchange field Bex, cr is the vector of Pauli matri- 
ces, and V describes the effect of potential scattering and 
spin-orbit scattering. In a ferromagnet, Ez is comparable 
in magnitude to the chemical potential /x. 

The free energy F is written as a contribution Fq in 
the absence of the potential V, which is isotropic, and a 
shift A_F. The shift AF is calculated as 

AF = r^trGo(c^„,k)Vk,k + ^T ^ trGo(c^„,ki) 

k,n ki ,k2 ,n 

X Vki,k2Go(t^„,k2)Vk2,ki + • ■ • , (3) 

where ujn — T:T{2n+ 1), n integer, is the fermionic Mat- 
subara frequency, Go{uJm^) is 2 x 2 matrix Green func- 



tion corresponding to Tio , and "tr " indicates a trace over 
the spin degrees of freedom. 

The Fourier transform yki.k2 of the potential V is 

Vki,k2 = Wk2-ki - ivll_^^(T ■ (k2 X ki), (4) 

where v and are the impurity potential and the spin- 
orbit potential, respectively. For the calculation of the 
correlator C (9) we need to perform the ensemble average 
over the disorder potentials v and v^°. We take these 
to be Gaussian random variables, as is appropriate for a 
small concentration of impurities, 



(5) 



The disorder strengths and are related to the 
golden-rule scattering times r and Tso for elastic scat- 
tering and spin-orbit scattering, respectively. We as- 
sume that disorder is weak and that scattering processes 
are predominantly spin-conserving, i.e., fJ,,Ez ^ H/t ^ 
h/TsQ. No assumption is made on the relative magnitude 
of the spin-orbit scattering rate H/tso and the mean spac- 
ing 5 between energy levels. For ease of presentation and 
for technical simplicity, the calculations below are for the 
case of a 'half metal', i?z ^ for which all electrons in 
the conduction band have the same spin direction. For 
the correlation function of the disorder potential we take 
the simple form 



= we{2kF ~ q), 



'e{2kF-q), (6) 



where kp is the Fermi wavevector of the half metal and 
0{x) = 1 if a; > and otherwise. In this simple model, 
the spin-orbit time for electrons at the Fermi level is given 
by Tso = QTrh^ /2w^°kpm, where hkp is the Fermi momen- 
tum. Although there will be quantitative differences, our 
order-of-magnitude estimates continue to hold for a true 
itinerant ferromagnet in which both minority and ma- 
jority electrons are present and if the disorder correlator 
differs from the simple model ©. 

Upon performing the disorder average for the correla- 
tor one finds that the leading contribution is from a 
series of ring diagrams, as shown in Fig. ^ Since there 
is no spin-orbit coupling for forward scattering, the first 
order ring diagram does not contribute to the mesoscopic 
anisotropy, whereas the higher order diagrams do. 

The second order contribution C2 is identical to what 
one would find in standard second order perturbation 
theory in V, except for the replacement of Gq by the im- 
purity averaged Green function. This replacement gives 
rise to a smearing of momenta by an amount of order of 
the inverse mean free path l/l, which is irrelevant for the 
final result. Hence, we find 
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C = 




FIG. 1: Diagrammatic representation of C{6). The dotted 
lines represent the impurity average. The solid lines are im- 
purity averaged single-electron Green functions. 



plus terms that do not depend on 9. The summations 
over q, ki, and k2 can be done analytically, with the 
result (again omitting terms that do not depend on 0) 



-Nc 



hsind 

27rTsn 



(8) 



where c = 3(5005 + 247r2)/25025 w 0.63 is a numerical 
constant and N is the number of electrons contributing 
to the magnetic moment. Below we show that the contri- 
bution to C (9) of the higher order diagrams is paramet- 
rically smaller than the second-order contribution C2 (9) . 
Hence, Eq. ((SJ contains the dominant contribution to the 
correlation function C{9). Higher order contributions are 
important only for the dependence of C (9) on an applied 
magnetic field or temperature However, such de- 

pendencies are weak and will not be considered here. 

It is important to note that the entire Fermi sea con- 
tributes to the magnetic anisotropy fluctuations. This is 
similar to, e.g., (unscreened) spontaneous dipole moment 
of a metal grain 13], but different from most other meso- 
scopic fluctuation effects, such as persistent currents and 
conductance fluctuations [l^ . which are determined by 
the properties of electronic states around the Fermi level 
only. As a consequence, C (9) is not universal; it depends 
on the details of the impurity potential or on the elec- 
tronic structure of the ferromagnetic particle. Indeed, in 
the general case of a ferromagnet with both majority and 
minority electrons, or for a more general disorder correla- 
tor than the simple model © , the value of the numerical 
constant c will be different, although the order of mag- 
nitude of C (9) and the angular dependence will be those 
of Eq. ||SJ) above. 

We now show that the contribution from the higher 
order diagrams to the correlator C{9) is small. Hereto 
we again consider the case of a half metal, Ez ^ /i. In 
order to remove the prefactor 1/n from the nth order 
diagram, n = 3,4, . . ., we calculate dC{9)/d9 instead of 
C(9). Because Tso ^ r, angular averages for the impu- 
rity lines and the single-electron Green functions can be 
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performed separately for the higher-order diagrams. Ne- 
glecting the g-dependence of the disorder correlator in 
Eq. lO, we then find 



q,p>0 ± 



de 



(9) 



where fip = 2nkTp is the bosonic Matsubara frequency 
and 



X±(f]p,q,0) 



Tso — t(1 ± C0S( 



Here the brackets (. . .)k denote an angular average over 
the direction of k. For small frequencies fip and wavevec- 
tors q one can approximate K as 



K±{np, q, 61) = 1 - VLpT - q^l^/Z 



(l±cos6'). (10) 



When either fip or q are large, one can approximate 



X±(r2p,q,( 



Tso — ^(1 ± cosO) 
rsomax(r2pT, 2ql / n) 



(11) 



The summation over p and q in the region in which both 
flpT < 1 and ^ 1 can be done using the approximation 
(|10|l for K. The zero mode (q = 0) gives a weakly singu- 
lar contribution to dC/d9, which is small in comparison 
to 02(0), whereas the nonzero modes give a non-singular 
contribution ~ V0/1^t^^ for small 9. The remaining part 
of the summation over p and q is done with the help 
of Eq. and gives a contribution to dC/d9 that is a 
factor ~ l/kpl ^ 1 smaller than C2. 

The above calculation fixes the magnitude of the meso- 
scopic anisotropy, but it does not fix the form of the 
anisotropy. In fact, both easy plane and easy axis 
anisotropics with a random orientation of the plane or 
the axis give the same angular dependence C{9) cx sin^ 9 
found from second order perturbation theory. The full 
angular dependence of F(m) compatible with C{9) cx 
sin^ 9 is one described by an anisotropy tensor _ff , 



a, [3=1 



(12) 



Our previous calculation of the free energy correlator has 
shown that anisotropy is dominated by second order per- 
turbation theory. Using second order perturbation the- 
ory to calculate the full tensor H for a half metal, we find 



-(kx q)a(kx q)^. (13) 



A similar but slightly more complicated expression holds 
for the general case. 



From a statistical point of view, the matrix structure 
of H is the same as that of a matrix Haf3 = '^ja'^jpi 
where the Uj are statistically independent vectors. The 
elements of such a matrix have Gaussian distributions, 
the variances of the diagonal elements being twice as 
large as the variance cr^ of the off-diagonal elements. 
Comparing the correlator C{9) — —2a^ sin^ 9 one obtains 
from a Gaussian distribution of H with the microscopic 
calculation performed above, one finds cr^ ~ Nfi^/T^^. 

Changing to principal axes, H can be brought to diag- 
onal form. Denoting the eigenvalues of H hy hi < h2 < 
h^, with corresponding eigenvectors vi, V2, and V3, one 
can rewrite Eq. 1)12(1 as 

F{m.) = -ha{m ■ vif + hp{m ■ vjf + const, (14) 

where ha ~ — hi and hp — h^ — /12 are anisotropy ener- 
gies for easy axis and easy plane anisotropy, respectively. 
The eigenvalues and eigenvectors of a real symmetric ma- 
trix with random Gaussian elements have a distribution 
that is well known from random matrix theory : The 
orientation of the eigenvectors Vi, V2, and V3 is fully 
random, whereas the distribution of the three eigenval- 
ues hi < h2 < is 



P{hi,h2,h3) (xY[{h, - h,)Y[, 



(15) 



This implies that the distribution of ha and hp is 



p^l^^ ) ^ ^hahpjha + fep)V6^ ^-(;»;+/,„/,,+/»g)/6^^ ^ 

(16) 

In particular, Eq. (|16|l predicts that the ratio r = ha/hp 
of easy plane to easy axis anisotropics has the distribu- 
tion 



P(r) 



27r(l 



8(r2 



1)5/2 



(17) 



This result is universal and does not depend on the details 
of the impurity potential or the electronic structure that 
affect the magnitude of the mesoscopic anisotropy. 

It is instructive to compare the magnitude and the an- 
gular dependence of the mesoscopic anisotropy to that of 
magnetocrystalline anisotropy. For magnetocrystalline 
anisotropy, the symmetry of the lattice often implies 
a more complicated angular dependence than that of 
Eq. ((14(1 1]. The magnitude of the magnetocrystalline 
anisotropy scales cx hN/Tgo, but the numerical prefac- 
tor may be small for polycrystalline samples, disordered 
ferromagnetic alloys, or ferromagnets with a close-to- 
parabolic dispersion relation. A crude comparison of 
magnetocrystalline and mesoscopic anisotropy for tran- 
sition metal magnetic nanoparticles can be made using 
the following estimates: Tso can be taken comparable to 
that in Cu [ig, Tso ~ 10~^^s, whereas the magnetocrys- 
talline anisotropy energy is 0.04 meV/atom for Co and 
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significantly smaller for Fe and Ni 0- With these esti- 
mates, the mesoscopic contribution to the anisotropy en- 
ergy should become dominant for Co nanoparticles with 
< 10^ atoms. For Fe and Ni mesoscopic anisotropy fluc- 
tuations start dominating already at larger particle sizes. 

Sofar we have dealt with the anisotropy energy F(m) 
at a fixed number of electrons. This energy determines 
the size of the energy barrier for magnetization rever- 
sal and the low-energy particle-hole and magnetic exci- 
tations. An important question for the experiments of 
Refs.i2,l3i|is how much the anisotropy constants ha and hp 
change upon addition of a single electron to the ferromag- 
netic particle. A theoretical estimate of 5ha and 5hp was 
given in Refs. 0, 0. Both theories assumed that there 
was a systematic source of uniaxial anisotropy, which is 
appropriate for the non-spherical nanoparticles used in 
the experiment; Addition of an electron caused a random 
change 5ha ~ h/T^o ^ ha/N 0. Here we briefly address 
the same question for a ferromagnetic nanoparticle with 
'mesoscopic' anisotropy only. As we'll show below, both 
5ha and 5hp remain of order h/Tso ^3|- However, since 
the 'mesoscopic' anisotropy is weaker, the relative magni- 
tude of the changes 6ha and 5hp is larger: 6h ~ hjN^I"^ . 

The change of anisotropy energy upon addition of a 
single electron in an energy level e follows from the de- 
pendence of £ on the magnetization direction m |^ 0| . 
The m-dependence of a single level e is well described 
by random matrix theory. For a half metal, there are no 
spin-flip processes, and one finds that e is an eigenvalue 
of the M X M hermitian matrix 



H(m) =S + i 



nh 



■ 3 



(18) 



Here 6 is the mean level spacing, 5 is a real symmetric 
matrix, and the Aa are real antisymmetric matrices, a = 
1, 2, 3. The elements of all matrices are chosen as random 
Gaussian variables, with the same variance for the off- 
diagonal elements. Diagonal elements of S have a double 
variance. The matrix size M has to be taken to infinity. 

In order to find how the addition of a single electron 
changes the anisotropy energy F(m) for m close to the 
direction mo at which _F(m) is minimal, we expand e(m) 
around mo. Statistical estimates of the first and second 
derivatives of e to m are relatively straightforward be- 
cause H{m.o) and _ff(m) — iJ(mo) may be considered sta- 
tistically independent for small angular deviations. The 
first and second derivatives of e to m are zero on av- 
erage. The mean square average of the first derivative 
is (|9e/9mp) = h6/2TTTso- Typically, the second deriva- 
tive is of order h/T^o- Using Eq. (|14|l for F(m), with 
ha, hp ~ hN-^/^ /tso, one finds that addition of an electron 
causes a shift of the preferred magnetization direction vxiq 
by an angle of order {TgoS/Nhy^^. The shifts Sha and 
dhp of the anisotropy constants for easy-axis and easy- 
plane anisotropy are ~ h/Tso, in agreement with Refs. 



0, 0. Note that the expansion of e around mo is valid 
up to an angle ^max ~ min[(rso(5/fi.)^/^, 1]. Both the an- 
gular shift of mo and the angular distance of different 
quantized levels are much smaller than ^max, so that the 
expansion used to obtain these estimates is justified. 

In conclusion, we have studied magnetic anisotropy in 
small ferromagnetic particles without shape or magne- 
tocrystalline anisotropy. The presence of non-magnetic 
impurities and spin-orbit scattering provides a 'meso- 
scopic' lower bound to the particle's anisotropy energy. 
This 'mesoscopic' magnetic anisotropy is caused by spin- 
orbit coupling, which connects the orbital anisotropy 
from non-magnetic impurity scattering to the spin degree 
of freedom. The strength of the 'mesoscopic' anisotropy 
is of order N^^'^H/tso, where N is the number of electrons 
contributing to the magnetic moment and Tgo is the spin- 
orbit time, the proportionality constant depending on the 
details of the electronic dispersion relation and the im- 
purity potential. The angular dependence of the 'meso- 
scopic' anisotropy is a random combination of easy-plane 
and easy-axis anisotropics. 
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